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Counting Lattice Triangulations 

Volker Kaibel and Giinter M. Ziegler 

Abstract 

We discuss the problem to count, or, more modestly, to estimate 
the number f(m,n ) of unimodular triangulations of the planar grid of 
size m x n. 

Among other tools, we employ recursions that allow one to compute 
the (huge) number of triangulations for small m and rather large n by 
dynamic programming; we show that this computation can be done in 
polynomial time if m is fixed, and present computational results from 
our implementation of this approach. 

We also present new upper and lower bounds for large m and n, 
and we report about results obtained from a computer simulation of 
the random walk that is generated by flips. 


1 Introduction 

An innocent little combinatorial counting problem asks for the number 
of triangulations of a finite grid of size m x n. That is, for m,n > 1 we 
define P m ,n '■= { 0 , 1 ,... , m} x { 0 , 1 ,..., n}, “the grid”. Equivalently, the point 
configuration consists of all points of the integer lattice 7? in the lattice 
rectangle conv(P m ri ) = [0,m] x [0,n] of area mn. Every triangulation of this 
rectangle point set that uses all the points in P m) „ has (m + l)(n + 1) = \P m , n \ 
vertices, 2 mn facets/triangles, and 3 mn + m + n edges, 2(m + n) of them on 
the boundary, the other 3 mn — m — n ones in the interior. All the triangles 
are minimal lattice triangles of area | (that is, of determinant 1), which are 
referred to as unimodular triangles. The grid triangulations that use all the 
points are thus called unimodular triangulations. The number of unimodular 
triangulations of the grid P m ^ n will be denoted by f(m,n). 

As an example, our first figure shows one unimodular triangulation of the 
5x6 grid: 



m = 5 
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To get started, one notes that of course f(m,n ) = one discovers with 

pleasure that f(l,n) = ( 2 "), and one works out by hand with a bit of pain that 
/(2,2) = 64 and /(2,3) = 852. Furthermore, one observes that the special 
triangulations that decompose into m vertical strips of width 1 yield the lower 
bound 

f2n\ m 

f(m, n) > /(l, n) m = (1.1) 


so for larger m and n the numbers /(m,n) get huge very soon; for example, 
the bound yields /(5, 6) > 6 ■ 10 14 . 


It is equally interesting to study/enmnerate more general types of triangu¬ 
lations, such as triangulations of finite point sets that do not necessarily use all 
the points, triangulations of general convex or non-convex lattice polygons, tri¬ 
angulations of general position point sets, triangulations of higher-dimensional 
point sets, etc. None of these will appear here, but we refer interested readers 
to Lee [|22| and De Loera, Rambau & Santos Jl2| . 


Lattice triangulations are basic combinatorial objects, and they are fun¬ 
damental discrete geometric structures; so it is no surprise that they appear 
in various computational geometry contexts. However, lattice triangulations 
have also been studied intensively from different algebraic geometry angles. 
So, triangulations of a convex lattice polygon 


provide the data for Viro’s |JT| famous construction method of plane 
algebraic curves with prescribed combinatorics and topology, related to 
Hilbert’s sixteenth problem; 

appear in GeLfand-Kapranov-Zelevinsky’s |[T7| theory of discriminants, 
where “regular” triangulations are in bijection with the vertices of the 
secondary polytope of the point configuration, and 
model torus-equivariant crepant resolution of singularities for toric three¬ 
folds, where “regular” triangulations correspond to projective desingu- 
larizations; see e. g., Kernpf et al. [fLL], Chap. 3] and Dais [|T0| . 


The last two points pose the problem of counting or estimating the num¬ 
ber / reg (m, n) of regular triangulations of P m ,m that is, of triangulations of 
conv(P m n ) whose triangles are the domains of linearity for a piecewise lin¬ 
ear convex function (lifting function ); see also, for example, Sturmfels 0, 
Chap. 8], Ziegler |32|, Lect. 5]. Motivated by the toric variety considerations, 
one would like to know whether/that for large m and n most triangulations 
are non-regular. There is no proof, but a lot of evidence that this should be 
true (see, e.g., Table |2|). 

In this context, we must admit that despite the effort that has been put 
into studying this question (see e. g. Hastings [ 0 . Chap. 2]), it has not become 
clear what non-regular triangulations really “look like.” The “mother of all 
examples” is the whirlpool triangulation of the 3x3 grid: 
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Indeed, this unimodular triangulation is irregular, since a lifting function h 
would have to satisfy h(pi ) + h(qi-i) < h(pi- i) + h(qt) for i = 0,..., 3 (all 
indices taken modulo 4), implying the contradiction h(q 0 ) + ••■ + h(q 3 ) < 
%o) + • • ■ + h(q 3 ). 

For m = n = 3, except for three symmetric copies this is the only example of 
a non-regular triangulation: We have /(3, 3)—/ reg (3, 3) = 4. This may lead one 
to the conjecture that some kind of “generalized whirlpools” are responsible 
for non-regularity. However, the pictures of “non-regular triangulations” that 
we present below do not support this intuition. 


The plan for this paper is as follows: In Section [2| we face the challenge to 
count explicitly, trying to cope with the “combinatorial explosion.” For this, 
we present a simple dynamic programming technique by which we get sur¬ 
prisingly far, and which on the specific problem of grid triangulations outruns 
the much more sophisticated general techniques such as Avis & Fukuda’s || 
reverse search, Aichholzer’s jl] path of a triangulation, or the oriented matroid 
technique of Rarnbau p7 |. 

In Section [|we estimate f(m,n) for large m and n, trying to narrow the 
bounds for the asymptotics. It is interesting to compare with the situation for 
N = (m + l)(n+1) points in the plane in general position, where the currently 
best available upper bound seems to be Santos & Seidel’s |28[ estimate that 
there are 0(59^) triangulations. However, in our problem the N points are not 
at all in general position — and the upper bounds that we have to offer are 
much better: We report on a neat 0(2 3iV ) upper bound by Anclin ||, which 
substantially improves on a previous 0(2 4iV ) upper bound by Orevkov 
Based on explicit enumeration results from Section ^]. we get a lower bound of 
22.055 mn w } ien both m and n get large; note that ( |1 . 1| ) yields already a lower 
bound 2 d-°(i)) 2 mn_ 


Finally, in Section [| we sample lattice triangulations for large parameters 
m and n, and thus try to understand what typical lattice triangulations, as 
well as typical regular lattice triangulations, “look like.” While we have some 
pictures to offer, proofs seem harder to come by. Indeed, the pictures display 
some long-range order; while this may make lattice triangulations interesting 
as a statistical physics model, it generates serious obstacles for any proof that 
the obvious Markov chain is rapidly mixing, and thus to application of the (by 
now) standard theory ||. 
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2 Explicit Values 


There are several methods available to generate or count all triangulations 
of a finite set of points in R 2 . An approach that works for point sets in arbitrary 
dimensions is implemented in the software package TOP COM by Rarnbau [27j 
(see also Pfeifle and Rarnbau |^ 6 |). It enumerates all triangulations in a purely 


combinatorial manner after the chirotope of the point set (oriented matroid 
data) has been computed. 

The reverse search algorithm proposed by Avis and Fukuda || is a rather 
general enumeration scheme that can be specialized to triangulations of point 
configurations in R 2 (see also Bespamyatnikh 0). Since it was used to ob¬ 
tain some of the results reported in Section [|, and because it is based upon 
some structural properties that are relevant for our treatment later, we briefly 
describe the method here. 

Let T be a fine triangulation of a point set S' C R 2 ; i.e., T is a set of 
triangles, for which the set of vertices equals S', such that the union of all 
triangles is the convex hull of S', and any two triangles intersect in a common 
(possibly empty) face. The unimodular triangulations of P m ,n are precisely its 
fine triangulations. An edge of some triangle in T is flippable if it is contained 
in two triangles of T whose union is a strictly convex quadrangle. Replacing 
these two triangles by the two triangles into which the other diagonal cuts 
that quadrangle ( flipping the edge) yields another fine triangulation of S'. The 
graph on the fine triangulations of S' defined via flipping is the flip graph of S'. 

Let us fix an arbitrary ordering of the points in S', inducing via lexicograph¬ 
ical ordering a total order on the set of pairs of points, and thus, again via 
lexicographical ordering, a total order on the fine triangulations of S', which 
are identified with their sets of edges here. There is a distinguished fine trian¬ 
gulation To of S' with respect to that ordering, namely the smallest Delaunay 
triangulation (i.e., a triangulation characterized by the condition that for every 
triangle the circumcircle contains no point from S' in its interior). 

Furthermore, for each fine triangulation T 7 ^ To, there is a distinguished 
flippable edge (computable in 0(|Sj) steps) such that, starting from any fine 
triangulation of S', iterated flipping of the respective distinguished edges even¬ 
tually yields Tq. This algorithmically defines a spanning tree in the flip graph 
of S', rooted at %; in particular, the flip graph of a two dimensional finite point 
set is connected. 

The basic idea of the reverse search method is to traverse that spanning 
tree from its root To. At each iteration one chooses a leaf T of the current 
partial tree, and determines those among the triangulations adjacent to T 
on whose path to T 0 in the spanning tree T lies. Properly implemented, the 
reverse search algorithm generates all fine triangulations in O(|S'! • f(S )) steps, 
where /(S') is the number of fine triangulations of S'; see ||. 

Via the “secondary polytope” of a finite point set S' C I 


16 


one can 


design a variant of the reverse search algorithm that generates all regular 
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triangulations of S in 0(| S\ ■ F Teg (S )) steps, where F Teg (S) is the number of 
such triangulations. It is, unclear, however, if one can also generate all regular 
fine triangulations of S in a number of steps that is bounded by a polynomial 
in the number of such triangulations. 

If one is interested in the number of triangulations of a two-dimensional set 
of points rather than in the explicit generation of all of them, then the path of 
a triangulation method due to Aichholzer [[!} is more efficient than the reverse 
search algorithm. 

For counting the unimodular (fine) triangulations of the very special point 
sets Pm,m however, different methods are much more efficient. These are de¬ 
scribed in the following sections. 

2.1 Narrow Strips 

Strips of width m = 1. For any lattice trapezoid of width 1, whose paral¬ 
lel vertical sides have lengths a and b, the number of unimodular triangulations 
isgi(a,b) = ( a ^ b ) = ( a ^ b ); indeed, a bijection between these triangulations and 
the a-subsets of {1,..., a + b} is established by top-down numbering the tri¬ 
angles of a triangulation T by 1,..., a + 6, and mapping T to the a-set of all 
numbers of triangles whose vertical edges are on the left. In particular, we 
have 



( 2 . 1 ) 


Strips of width m — 2. For /(2,n) we have no explicit formula, and we 
cannot evaluate the asymptotics precisely, but we have a “quadratic” recursion 
that can be evaluated efficiently: For this we enumerate the triangulations 
according to the highest “width 2 diagonal,” which (if it exists) decomposes 
the rectangle into two width 1 strips, a single triangle, and a trapezoid of 
width 2 (see the left-hand figure below). Let g 2 (A, B) denote the number 
of triangulations of a trapezoid of width 2 with vertical edges of lengths A 
and B and horizontal base line, as in the right-hand figure below — where 
A + B = 1 mod 2 implies that the midpoint of the diagonal is not a lattice 
point, and where we may assume A < B by symmetry: 


trapezoid shape 
(A,B): 




A 


A 
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Thus we get 
/(2,n) = 


2 n 


n 


2 £ 9 ,(4B)( 2n ; 


3A+B+1 


0<A<B<n 
A-\-B =1 mod 2 


A 


2n 


n 


A+3-B+1 

2 

-5 


The binomial coefficients in this recursion correspond to triangulations of 
width 1 strips. So they could be rewritten in terms of gi, as gi(n—A, n— A+ ^ +1 ) 
resp. gi(n — B, n— A+ ^ +1 ). A similar remark applies to the binomial coefficients 
that appear in the following. 

For (A B) = g 2 (B, A) we also get a recursion by considering the highest 
diagonal of width 2: 

/ 3A-t-B-l \ / A+3B-1 ' 

g 2 {A,B) = 


2 

A 


2 

B 


E 


92(a,b) 


0 <a<A, 0 <b<B 
a-\-b= 1 mod 2, a+6<A+S 


3A+B—3a—fc 
2 

A — a 


A+3£-a-3fe 

2 

B-b 


A 


B 


A 


B 


Here the parameters A, B, a, b may be interpreted as the ^-coordinates of cer¬ 
tain lattice points. The shaded parts of the figure consist of strips of width 1, 
whose triangulations are counted by the binomial coefficients gi(-, •). 


Strips of width m = 3. For /(m, 3) we have a recursion of order 4; 
it relies on the observation that if we screen the middle strip from the top 
for diagonals of width at least 2, then the first diagonal to find will be of 
width exactly 2, since any width 3 diagonal is flippable, and contained in a 
parallelogram that is bounded by two width 2 diagonals. The corresponding 
decomposition of our rectangle is indicated in the left figure below: 
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Therefore, we obtain 

/( 3,n) = ( 2 ™) +2 h(A,B,n, 

' ' 0 <A,B<n 

A-\-B =1 mod 2 


2 n 


3A+B+1\ 


n 


2 

A 




2n 


n 


A+3B+1 

2 

5 


where h(7l, B , C, D) counts the number of triangulations in a “hook” shape as 
given in the right drawing in the figure above, which depends on four param¬ 
eters 0 < A, B,C, D < n, with A + B = 1 mod 2 and B < C. For the number 
of triangulations of such a hook shape we get a recursion 


h{A, B , C, D ) = 


3A+B-1\ /A+3B-1 


2 

A 


2 

5 


C + D 
C 

/ 3A+B—3a—b 


+ 


E h(a,b,C,D)i 2 

4 n<rh<TR \ 


0<a<A, 0<b<B 
a-\-b= 1 mod 2, a+b<A-|-S 

E 

0 <a<D, 0 <b<A+B=l 
a-\-b=l mod 2, <_B 


A+3B-a-3b 

2 

B-b 


h(a , b, 


A+B—l 


,A) 


D + C 


3a+fc+l 


D — a 


A+3B~a~3b 


- 1 


B 


ci-(-6-1-1 

2 




3.B-A-1 A+B-l A+B-l 




C + D 


5-B—A—1 


C-B 


if D > > 0. 


The four terms in this recursion correspond to the four cases depicted in the 
figure below, where the fourth case — of a long diagonal of width 3 — occurs 
only in the case where the second endpoint of the diagonal, which may be 
computed to have ^-coordinate , comes to lie within the hook. 



2.2 Strips of (fixed) width m. 

We now describe a recursive strategy for the enumeration of unimodular 
triangulations of grids of arbitrary size. The method is applicable for triangu¬ 
lations of general finite point sets — but it is effective only in the special case 
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where the points lie on a small family of parallel (vertical, say) lines; in our 
case this is the situation of small (fixed) m and variable n. The key observation 
is that any triangulation may be dismantled by removing triangles from the 
upper boundary, while maintaining a lattice triangulation of a y-convex lattice 
polygon. Since for fixed m the number of such polygons in P m , n is bounded by 
a polynomial in n, this yields an efficient dynamical programming algorithm 
in this case. 

Let A 1: A 2 C R 2 be two triangles whose intersection A x n A 2 is a common 
(empty, zero-, or one-dimensional) face of both A x and A 2 . We say that 
A 2 lies above A x (A x -< A 2 ) if there are two points (x, yi) G A x \ A 2 and 
(x, y 2 ) G A 2 \ A x on a vertical line with y x < y 2 . For example, in our figure 
the shaded triangle lies above the other one; the other two pairs of triangles 
are incomparable. 



Due to the convexity of the triangles and the intersection condition imposed 
on them, this is a well-defined asymmetric and irreflexive relation. 

Lemma 2.1 There is no sequence A 0 ,..., A t -i C R 2 of triangles (such that 
the intersection of any two among them is a face of both) satisfying 


Ao -< A x -<...-< At-i -< Ao . (2-2) 

Proof Suppose that A 0 ,..., A f _ x C l 2 is a minimal cyclic sequence of tri¬ 
angles (such that the intersection of any two among them is a face of both), 
i.e., it satisfies (|2.2|) (it is cyclic ), but no subsequence of A 0 ,..., A t _ x C R 2 is 
cyclic. The orthogonal projections x(Aj) to the x-axis have the following three 
properties (where all indices are taken modulo t ): 

(a) x(Ai) Dx(A i+1 ) ^ 0, 

(b) x(A i) <£ x(A i _i),x(A i+ i), and 

(c) x(Aj_i) n x(A m ) = 0, 

where (a) follows immediately from the definition of -< and (b) as well as (c) 
are due to the minimality of the cycle. 

But (a), (b), and (c) together imply that the intervals x(A 0 ),... , x(Af_ x ) 
“either run left-to-right or right-to-left”; in particular, we have x(A 0 )nx(Aj) = 
0 for i G {2, — 1}, contradicting A-< A 0 . □ 


Of course, both in the definition of -< as well as in Lemma |2T 
replace “triangle” by “compact convex set.” 


one can 








Counting Lattice Triangulations 


9 


The relation -< was defined with respect to parallel projection here. If one 
defines an ordering with respect to central rather than to parallel projection, 
then the analog to Lemma 0 (for arbitrary centers of projection) does not 
hold. For Delaunay triangulations, however, there is an analog to Lemma 
(See De Floriani et al. 
dimensions.) 


lfil for dimension 2, and Edelsbrminer [[L3(] for arbitrary 


Due to Lemma 2.1 

p2 


the relation -< induces a partial order on the set of 
triangles in M 2 , which we will also denote by -<. 

A sequence 7),..., P 2 mn of sets of triangles will be called an admissible 
sequence for P m ^ n if Tj is a unimodular triangulation of P m ,n, and if, for each 
i = 2,3,..., 2 nm, we have % = %-\ \ {A} for some -<-maximal triangle A 
in A subset S' C M 2 is called an admissible shape (of P m , n ) if it can 

be obtained as a union S = |J Agr . ^ f° r an admissible sequence 7j,..., T 2rnn 
and some i G {l,...,2mn}. Every admissible shape is y-convex (i.e., its 
intersection with any vertical line is connected). It is determined by its upper 
boundary segments , i.e., the sequence of line segments [/T)^! 1 )],... ; 
with lx' 1 = 0, lx } = n, and for j G {2,... ,t}, such that, for each 

point p in the relative interior of any of the segments, p + (0, e) S holds for 
all e > 0. 

Let S be an admissible shape. We denote by 7 max (S') the set of all 
maximal unimodular triangles in S, that is, the finite set of all unimodular 
triangles that could be -<-maximal in some unimodular triangulation of S. 
For example, the figure below indicates the 12 -<-maximal triangles of the 
shaded admissible shape. 



Any admissible shape S' that arises from S by removing some triangles con¬ 
tained in 7j nax (S') is called an admissible subshape of S. These triangles have 
disjoint interiors, and they are uniquely determined for each admissible sub¬ 
shape S' of S (compare the proof of Lemma |2.5| ). Their number is denoted by 
#(S',S). 

Since every unimodular triangulation of S contains at least one triangle 
from 7j nax (S'), we obtain the following inclusion-exclusion formula for the num¬ 
bers f(S) of unimodular triangulations of admissible shapes S. 
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Lemma 2.2 Every admissible shape S has 


f(s ) = 5J-I) #(s, ' s) -7(S') 


(2.3) 


S' 


unimodular triangulations, where the sum is taken over all admissible proper 
subshapes S' of S. 

Lemma ^[2] allows us to compute f(m, n ) = /([0, m] x [0, n]) via a dynamic 
programming approach: Determine the numbers /(S') via (|2.3|) in some order 
such that every admissible shape appears after all its admissible subshapes. 
In order to analyze the running time of such an algorithm, we first need to 
estimate the number of admissible shapes. 

Lemma 2.3 Let [/d) jr ( 1 )] ) ... ; j £/ ie sequence of upper boundary seg¬ 

ments of an admissible shape S. Then 



holds for each 2 < j < t. 

Proof This follows by induction on the number of triangles removed in order 
to obtain S : Every vertical edge of a triangle in a unimodular triangulation 
has length one, and thus removing a -X-maximal triangle from an admissible 


shape never creates a vertical boundary part of height more than 1. 


□ 


Lemma |2.3| implies the following bound on the number of admissible shapes. 

Lemma 2.4 There are at most (3n + 2) m ~ 1 (n +1) 2 admissible shapes of P mn . 

Proof The upper boundary of an admissible shape has n +1 possible start and 
n + 1 possible end points. At every interior x-coordinate (x G — 1}) 

either a segment of the upper boundary ends and a new one starts (3(n +1) — 2 
possibilities, by Lemma [2~3]) or a segment “passes through” (one possibility). 
□ 

The second important quantity for the analysis of the running time of the 
dynamic programming algorithm proposed above is the maximal number of 
summands that may occur in (|2.3|) . 

Lemma 2.5 Every admissible shape S of P m ,n has at most 



admissible subshapes. 
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Proof Let (B 1 ,..., B t ) be the sequence (from left to right) of upper boundary 
segments of S. Each triangle in T max (S) contains at least one of the edges 

Let B E {Bi, ..., B t } be a segment of the upper boundary of S. There are 
at most two triangles in T ma . x (S) containing B and one of its adjacent segments. 
Each other triangle in T max (S) that contains B must have its third vertex v 
in B + R • (0, —1) on the line that is parallel to B at distance ^(-B) -1 , where 
£ 2 (-£>) denotes the Euclidean length of B (because the triangle has area |). 
Since v must be integral and there is no integral point in the relative interior 
of B , there are at most two possibilities for v. Let us call one of them the first 
triangle below B, and the other one, the second triangle below B (if they exist). 

For every admissible subshape T of S, define a word w(T) E (0, a, ft, '•/}* 
by replacing Bi by 

‘cd if S' \ T contains the first triangle below B,, 

‘ft’ if S' \ T contains the second triangle below B,, 

‘ 7 ’ if S' \ T contains the triangle formed by Bi and Bj_i, and 
‘O’ otherwise. 


Clearly, every ‘ 7 ’ in w(T) has a ‘O’ as its left neighbor; furthermore, w(-) is 
an injective mapping. Therefore, the number of admissible subshapes of S' is 
bounded from above by the function (p(m) defined recursively via 

79 ( 0 ) = 1 , 79 ( 1 ) = 3 , (p(m) = 3 ip(m — 1) + tp(m — 2) (m> 2). 


Using standard techniques (see, e. g., |29], Thm. 4.1.1]), one derives from this 
recursion that ip(m) < (!+LH) m f or m > 1 □ 


Lemmas 2A and |2.5| imply that (for every m) the dynamic programming 
algorithm needs at most 

3.31 m (3n + 2) m-1 (n + l) 2 < 10 m (n + l) m+1 


arithmetic operations. The actual running time of an implementation heavily 
depends on the data structures for storing the admissible shapes and on the 
way by which one determines the admissible subshapes in that data structure. 
Therefore, here we include only the following rough statement. 


Theorem 2.6 For every fixed m, the function f(m,n ) can be computed in 
time bounded by a polynomial in n. 


In our implementation, we organize the admissible shapes as the leaves 
of a tree whose nodes are the prefixes of the sequences of upper boundaries 
segments of admissible shapes. The data structure allows quite efficient access 
to the admissible subshapes while not wasting too much memory. Nevertheless, 
the bottleneck in the computations is always memory. It is crucial to use an 
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ordering of the admissible shapes in which for each shape S the shapes of which 
it is an admissible subshape come as soon as possible after it; this allows one to 
keep in memory only a subset of the admissible shapes at each point of time. 

Furthermore, in ( |2.3| ) there is no need to sum over all admissible sub¬ 
shapes S'. It suffices to consider those S' that arise from removing triangles 
hi Tmax(S') that are contained in {(x, y) G R 2 : x > c}, where c is the maxi¬ 
mal ^-coordinate where ly J> = ry -1 ^ + 1 for some j and the upper boundary 
segments [/( 1 ) ) , [l^\r^] of S (if that maximum exists). 

Of course, when computing the number of unimodular triangulations of 
Pm,n by our method, we obtain as a byproduct the number of unimodular 
triangulations of several interesting polygons inside P m ,n , including P m> i, ..., 
Pm,n— 1 • 

The algorithm described above cannot only be used to calculate the number 
of unimodular triangulations of any admissible shape S in P m , n ; it can also be 
extended to produce a uniformly distributed random triangulation within the 
same asymptotic running time, and thus, in polynomial time (depending on n) 
for fixed m. For this, one just determines the numbers of triangulations of those 
admissible subshapes of S that arise from removing one single triangle; with 
respect to the corresponding random distribution one then chooses one of the 
-^-maximal triangles in S at random, and proceeds with the subshape obtained 
by removing it. 


2.3 Explicit values 


We have implemented the algorithms described in Subsection |2.1| in C++, 
using the grap library |I7J for exact arithmetic, with the interface to it provided 
by the polymake system |15 . 

Results obtained by our code are compiled in the Appendix (see Tables ||, 
USB and |7|). The number of unimodular triangulations of P m)Tl asymptoti¬ 
cally grows exponentially with mn (see also Section ^|). Therefore, it is more 
convenient to view the function /(m, n) on a logarithmic scale, normalized 
by mn. 

Definition The capacity of the m x n grid is 

log 2 f{™, n ) 


c(m,n) : = 


mn 


The following Figure shows the capacity functions c(m, n) for m 6 {1,..., 6}. 
The largest capacity we found is 2.055792 (for m = 4 and n = 32, see Table |5|). 

To give an impression of the amount of (machine) resources required for 
the calculations: The run of the admissible shape algorithm for m = 6 and 
n = 7 needed about three gigabytes of memory. The grid P 6j 7 has 370252552 
admissible shapes. We generated them in lexicographical order with respect 
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to the pairs of starting heights and volumes. This way, never more than 15% 
of the admissible shapes had to be kept in memory simultaneously. Notice 
that more than 400 megabytes (of the 3 gigabytes in total) where needed 
just to store the numbers of triangulations of the admissible shapes in the 
memory. The CPU time used for the computation was about 25 hours. Our 
computations were performed on a SUN UltraAX MP machine equipped with 
four 448 MHz UltraSPARC-11 processors (of which we used only one) and 4 
gigabytes main memory. 

For very small parameters, Meyer [^4| has enumerated all unimodular tri¬ 
angulations by Avis and Fukuda’s reverse search method (sketched at the 
beginning of this section) and checked them for regularity. Table |l] shows the 
results. While for these small parameters irregular triangulations are quite 
rare, the picture changes drastically when m and n get larger, see Section [|. 


m x n 

# triangulations 

# irregular 

fraction 

3x3 

46456 

4 

.000086 

3x4 

2822648 

502 

.000178 

3x5 

182881520 

63528 

.000347 

4x4 

736983568 

1553020 

.002107 


Table 1: Number of regular triangulations of small grids. 


3 Bounds 

3.1 Patching 

Any two unimodular triangulations of P mm and P m ,n 2 can be patched to 
a unimodular triangulation of P m ,n 1 +n 2 - Thus we have the follwing super- 
multiplicativity relation, where / irreg (m,n) denotes the number of irregular 
unimodular triangulations of P mn . 
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Lemma 3.1 For m,ni,n 2 > 1 the following relations hold. 

(i) f{m,n 1 +n 2 ) > f(m,n 1 )f(m,n 2 ) 

(ii) / irreg (m, ni + n 2 ) > f(m, ni)/ irreg (m, n 2 ) 

With respect to regular triangulations, patching is dangerous, as demon¬ 
strated by the following example of a non-regular triangulation of P 4 4 com¬ 
posed of four regular triangulations of P 2)2 (suggested by Francisco Santos). 



However, a much more general theorem by Goodman and Pach [|18|] says that 
any two regular triangulations of two disjoint convex polytopes P 4 ,P 2 C 
can be extended to a regular triangulation of conv(P! U P 2 ) without additional 
vertices. Thus also for the regular case we get a (slightly weaker) supermulti- 
plicativity relation. 


Lemma 3.2 For m,ni,n 2 > 1 we have 


f res (m, n\ + n 2 + 1) > / reg (m, ni)/ reg (m, n 2 ) . 


The following figure illustrates the patching of Lemma 3.2 



For n 2 — 1 we will (Lemma 
Let us fix some notations first 


|) strengthen the inequality in Lemma |3T2 . 


For any function h : P t 


m,n 


we denote 


by H : 


• M 3 the function with H(x, y ) = (x, y, h(x, y)). The function h 
is called a lifting function of a triangulation T of P m>n , if T is the image of 
the set of “lower facets” of the 3-polytope con v{H(x,y) : (x,y) G P m . n } under 
orthogonal projection to the x, y-plane (deletion of the third coordinate). 

A function h is a lifting function of T if and only if h is convex and piece- 
wise linear, and its (maximal) domains of linearity are the triangles in T. In 
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particular, one may add to h any convex piecewise linear function whose do¬ 
mains of linearity are unions of triangles of T in order to obtain another lifting 
function for T. 

A triangulation is regular if and only if it has a lifting function. 

The following result shows that all unimodular triangulations of a strip 
of width 1 are regular; moreover one has a lot of freedom in choosing the 
respective lifting functions, which we will exploit below. 

Lemma 3.3 Let T be a unimodular triangulation of a lattice trapezoid with 
two parallel vertical or horizontal sides So and S± at distance one. Every 
piecewise linear function h 0 : So —> M that is strictly convex on So D Z 2 can 
be extended to a lifting function for T. 

Proof Let p 0 , ■ ■ ■ ,p r be the integral points on S 0 , and let {e i; i,..., e^-p)} be 
the edges of T connecting pi to S\. Let S 1 ^ (z) and Sf{i) be the closures of the 
two components of So \ {pi}- Then we may decompose the function h 0 as 

r k(i) 

ho(x) = • 

*=0 j =1 

where each h tJ is a convex function on S 0 , linear both on S'J'(z) and on Sf(i), 
and having its unique break-point at pt. 

Now we extend each hij to a convex, piecewise-lincar function defined on 
the entire trapezoid such that it has a break-line at the edge and is linear 
above and below this line. Then the sum of all h tJ is a lifting function for T. 

□ 


Proposition 3.4 For n > 1 the following relations hold. 

(i) /"*(l,n) = /(l,n) = CT) 

(W) r*(2,n) = /(2,n) 

Proof Part (i) follows immediately from Lemma |J(i) (and equation ( |2.1| )). 

Since patching two regular triangulations along a single edge preserves 
regularity, it suffices for the proof of part (ii) to show that every unimodular 
triangulation of shapes of one of the forms 
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is regular. But this can be derived from Lemma 3.3. One starts from an 


arbitrary prescribed strictly convex function on the middle column, and after 
the extensions to the two shaded vertical trapezoids of width one obtained 


from Lemma |3.3| , one adds a piecewise linear function that is constant on the 
left strip, linear on the right strip, and sufficiently large on the right column 
of vertices. □ 


Similarly, one proves the following strengthening of Lemma [L2] for 77-2 = 1, 
as announced earlier. 


Lemma 3.5 For m,n> 1, we have 


/ reg (m,n + 1) > / reg (m, n) • 


2 n 


n 


3.2 Limit capacities 

In the following, we will show that the capacities c(m,n) and c reg (m,n) 
(with f(rn,n) = 2 c ( m > n )" m and / reg (m, n) = 2 cies ( m ’ n l mn ) asymptotically behave 
well, which allows us to focus on their limits subsequently. Note that all 
capacities are bounded (see Theorem |3.9|). 


Proposition 3.6 Let m > 1. 

(i) The limit c m := lim c(m,n ) exists. 

n—>oo 

(ii) The limit c^ g := lim c reg (m,n) exists. 

71—KX) 


Proof Lemmas |3.1| (i) and |3.2| imply by Fekete’s lemma [23, Lemma 11.6] that 
lim and lim / reg (m, n — 1)« 

n —>00 n —>00 

exist. Therefore, 

r 1 1 f( \i r l °sKm,n) 

Inn —log j(m,n) n = inn -= inn c(m : n) 


n—>oo m 


rwoo mn 


and 

lim — --- log / reg (m, 71 «- 1)« = lim log f Tes (™,n - 1) = Um c re g ^^ 

n^oo m(n — 1) n^o o m[n — 1) n^oo 


exist as well. 


□ 


While the last proposition concerned the asymptotics of growing n for 
fixed m, the next result shows that also growing m and n simultaneously 
yields nice asymptotics. 
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Proposition 3.7 Let m > 1. 

(i) The limit c := lim c(m,m ) exists. It satisfies 

m—xx 

c = lim c m and c mo < c (m 0 G N) . 

ra—xx) 


(ii) The limit c reg \= lim c reg (m,m) exists. It satisfies 

ra—> oo 

c reg = lim c r J and c r ^ g < c reg (m 0 G N) . 

m—xx 0 


Proof From Lemma |3T| (i) one derives (for mo, no > 1) the inequality 

/(m,n) > (3.1) 

For integers p, q > 1 we define <3>(p, g) := 1 — p mod g . We have <f>(g, g) = 1 and 


lim 


XX) 


$(p, g) = 1 for all g G N. 


Equation ([LTJ) then implies 

c(m,n) > $(m, m 0 )d ) (n, n 0 )c(m 0 , n 0 ) , 

and. in particular, 

c(m 0 ,n) > <F(n,no)c(m 0 ,no) . 
Inequality ( |3.3| ) (together with lim <F(n, no) = 1) yields 


Inequality 


and therefore, 


c m > c(m,n 0 ) ■ 

implies (together with lim $(m,mo)$(m,no) = 1) 

ra—xx) 

lim inf c(m, m ) > c(m 0 , n 0 ) , 


(3.2) 

(3.3) 

(3.4) 


lim inf c(m, m) > c 


mo 


(3.5) 


Finally, we obtain the following chain of inequalities, which, together with 
( fi.5| ) proves part (i) of the proposition. The middle inequality is from (|3.5|), 
and the outer ones are due to (3.4). 

lim inf c m > lim inf c(m, m ) > lim sup c m > lim sup c(m, m ) 

ra—xx) ra—xx) m _ ) . 00 ra—>-oo 

Part (ii) is proved similarly, starting from Lemma f>.2| . □ 

Note that a similar proof yields 


c = lim c(am, j3n) and c reg = lim c Te& (am,f3n ) 

ra—xx) m—xx) 

for each pair a, /3 > 0. 

By Lemma |3.1| (ii), the corresponding “irregular limit capacities” do exist 
as well, and they are equal to c m (m > 3) and c, respectively. Therefore, we 
do not treat them explicitly. 
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3.3 Lower bounds 

Proposition 3.8 The following estimates hold. 

(i) Cl = = 2 

(ii) c 2 = c r 2 eg > 2.044 

(iii) c 3 > 2.051, c 4 > 2.055 

(iv) c m > 2.048 for m > 5 

(v) c m > c™* > 2 for m> 3 

Proof Part (i) follows from 

/(l,n) = / reg (l,n)= Im¬ 
parts (ii) and (iii) are results of the computer calculations reported in Section |. 
combined with Proposition |3.4| (ii). 

Lemma ^Tl](i) applied to the “transposed grids” implies 

\ m l / 777,2 , 

c(mi + m 2 , n) >-c mi, n) H-c(m 2 , n ) , (3.6) 

777,1 + 7?7 2 777i + 7?7, 2 



and thus 

777,1 777, 2 

Gni+m2 — ' C mi "h " C, 

777-1 + 777 2 777-1 + 777 2 

With (ii) and (iii), this yields (iv). 

Similarly, Lemma |3.5| leads to 


reg 

L 'm +1 


> 


777 

777 +1 


^ re g 


+ 


l 

777 +1 


• 2 . 


With (ii) this proves part (v). □ 

Equation (|+6|) implies that Ck n > c n for k > 2; for example, this implies 
that C\ > c 2 — but it is not obvious that c 3 > c 2 . Even stronger, one would 
assume that 

2 = ci < c 2 < c 3 < • • • < c , 

and 

2 = ci eg < c r 2 eg < c^ eg < • • • < c reg , 

but neither monotonicity is proved. 









Counting Lattice Triangulations 


19 


3.4 Upper bounds 

From general principles (see Ajtai et al. [[|]) one gets that the capacity for 
any configuration of N points in the plane is finite. In the general position 
case (no three points on a line) the currently best upper bound is o(2 59iV ), 
due to Santos & Seidel [^]. However, in the very “degenerate” case of lattice 
triangulations, there are far fewer triangulations: Orevkov |^5] obtained the 
bound f(m, n) < 4 3mri = 2 6mn . Very recently, this has been substantially 
improved by Anclin ||, as follows. 

Theorem 3.9 (Anclin [^]) For all m,n > 1, 

f(m, n) < 2 3mn-m-n . 


Proof Our sketch relies on the essential ideas of Anclin’s proof. 

The first, crucial observation is that the midpoints of the edges in any 
unimodular triangulation are exactly the half-integral, not integral points 
in coiw(P m>n ). (Clearly the midpoint of every edge is half-integral; the con¬ 
verse may be derived from Pick’s theorem, or from the fact that all unimodu¬ 
lar triangles are equivalent to conv{(0, 0), (1, 0), (0,1)}.) The number of these 
half-integral points in the interior of P m ,n is e = 3 mn — m — n. 

Now any triangulation is built as follows: The half-integral points are pro¬ 
cessed in an order that is given by a parallel sweep. (See de Berg et al. || 
Sect. 2.1] for a discussion and many further sweeping applications of this fun¬ 
damental technique.) 

Whenever a point is processed, we add to the partial triangulation a new 
edge with the given midpoint. The key claim is that at each such step, when 
a half-integral point v is processed, there are (at least one and) at most two 
possibilities for the new edge with midpoint v to be added. If this claim is 
true, then the number of triangulations is bounded by 2 e . 



To prove the claim, one can verify the following: Let [r/,1/] and \v" ,v"\ be two 
potential edges with midpoint v that could be added, where v' and v" are the 
endpoints below the sweep-line i\ let Q be the convex hull of the integral points 
in the triangle [v,v', v"]; then all the k > 2 vertices v' = vo,vi,, Vk-i, Vk = v" 
of Q are “visible” from v, and any one of the edges [vi, Vi] with midpoint v could 
potentially be added at this step. Furthermore, the midpoints |(uj_i + Vi) lie 
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below the sweep-line, and thus the edges are present in the current 

partial triangulation. 

v 2 v 0 



Now assume that [n 0 , w 0 ], [ui,hi], [u 2 ,h 2 ] are three adjacent edges with mid¬ 
point v that could be added when processing v. By central symmetry with 
respect to v we may assume that the midpoint w of [fi,h 2 ] lies below the 
sweep-line. But the triangle [vi, v , v 2 ] is then an empty triangle of area just 
like the triangle [ui,u,v 2 ]. From this we conclude that in the current partial 
triangulation, the edge [r>i, v 2 \ must be present — which creates a crossing with 
the potential edge [uo,ho], and thus a contradiction. □ 

The following upper bounds on the limit capacities follow immediately from 
Theorem |3T| . 


Corollary 3.10 For all m,n>l, the following inequalities hold: 

(i) c reg (m, n) < c(m,n) < 3 — ^ 

(ii) < c m < 3 — d- (in particular, cf e < c 2 < 2.5) 

(iii) c reg < c < 3 

As Anclin noted, his proof works much more generally: For any partial 
triangulation of a not necessarily simple or convex lattice polygon, the number 
of completions is at most 2 e , where e' is the number of edges that are to be 
added. 

4 Explicit Triangulations 

For small grids, one can enumerate all unimodular triangulations by the re¬ 
verse search algorithm sketched in Section ^]. For larger grids, it is desirable to 
obtain from the huge set of unimodular triangulations “random” ones. There 
are ways to produce them, however, in most cases the probability distribution 
from which they are chosen is unknown. 
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4.1 Generating Random Triangulations 

A standard way to compute complex random objects such as triangulations 
is to set up a random walk. In our case of unimodular triangulations of P m>n the 
method is described easily: First one determines any starting triangulation T 
of Pm,n- Then, the following operation is performed r times: Choose one 
of the (inner) edges of T uniformly at random; if this edge is flippable (see 
Section 2) , then with probability \ the current triangulation T is replaced by 
the one obtained from it by flipping that edge. 

As the flip graph of the triangulations is connected (see Section H), it 
follows from general principles that, with r tending to infinity, the probability 
distribution defined by the output of this algorithm converges to the uniform 
distribution (with respect to the “total variation distance”); see, e.g., Jerrum 
and Sinclair |2(J or Behrends [§]. Unfortunately, not much is known about the 
speed of convergence of the distribution. In particular, for general m and n 
it is not known whether there is a polynomial bound (in n + m + e _1 ) on the 
number r of steps needed to guarantee that the total variation distance between 
the produced and the uniform distribution is at most e (i.e., if the associated 
Markov chain is rapidly mixing). The only exception is the case m — 1: Here 
it follows from results of Felsner and Wernisch [jnj] that the Markov chain is 
indeed rapidly mixing. 

Despite this lack of knowledge on the distribution of the output of the 
random walk algorithm, one still can use it in order to produce examples of 
“interesting” triangulations. 


4.2 Empirical Results 


For each n € {10,... ,20}, Meyer |^4[ generated 1000 random unimodular 
triangulations of P nn by running the random walk described in Subsection [4T 
10 9 steps, recording every 10 6 th triangulation; see Table [2|. 

The results support the conjecture that, for n tending to infinity, a random 
unimodular triangulation of P n , n is irregular with probablity one. A second 
observation is that the (expected value) of the average length of an edge in a 
random triangulation seems to grow very slowly with n. 


4.3 Obstructions to Regularity 


All figures shown below have been produced by Meyer |23|, who imple¬ 
mented procedures for checking regularity by solving linear programs using 
the CPLEX 6.6.1 library. 

Proposition 373 shows that P^^n has only regular triangulations. The grid 
P 3>3 has precisely the following four (pairwise congruent) irregular unimodular 
triangulations. 
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m x n 

irregularity 

max. edge length 

av. edge length 

10 x 10 

.355 

5.538 

1.614 

11 x 11 

.435 

5.843 

1.630 

12 x 12 

.559 

6.118 

1.645 

13 x 13 

.696 

6.397 

1.659 

14 x 14 

.782 

6.650 

1.670 

15 x 15 

.875 

6.911 

1.681 

16 x 16 

.927 

7.151 

1.690 

17 x 17 

.965 

7.391 

1.700 

18 x 18 

.971 

7.618 

1.708 

19 x 19 

.992 

7.821 

1.713 

20 x 20 

.997 

8.060 

1.723 


Table 2: Results for random unimodular triangulations of large grids. The 
first column shows the (empirical) probability of irregularity, the second and 
third columns contain the (empirical) expected values of the maximal and the 
average edge length. 



When trying to understand the reasons for irregularity, it seems useful to 
consider (smallest) forbidden patterns for regular triangulations. Let T be a 
set of unimodular triangles of P m ^ n (such that any two of them intersect in 
a common face). We denote by S C P m ,n the set of all grid points covered 
by triangles in T. The set T is called regular if there is a height function 
h : S —> M such that for each triangle A G T, all /i-lifted points in S \A lie 
strictly above the affine hull of the h-lifting of A. A subset of T is called a 
minimal irregular configuration if it is not regular, but all its proper subsets 
are regular. Clearly, a regular unimodular triangulation of P m ,n cannot contain 
any minimal irregular configuration. 

Examples for minimal irregular configurations of P 33 are the following: 



Already for P 3 4 many other minimal irregular configurations occur; some 
are depicted here: 
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While these figures still have some similarities with the nice “whirlpools” 
ones for P 3 3 , the picture gets more and more complicated with growing grid 
sizes, as the following examples demonstrate: 



Viewing these figures, it seems unlikely that one can find any compact 
characterization of regularity for unimodular triangulations of P m: n in terms 
of forbidden substructures. 

We close our zoo of “explicit triangulations” with some pairs of triangula¬ 
tions, found by Meyer’s implementation of the random walk. In each of the 
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figures, the left triangulation is regular, but the right one is not, although it 
can be obtained from the left one by flipping just one edge (drawn bold in 
the upper left and in the lower right corner, respectively). For both irregular 
triangulations, a minimal irregular configuration contained in it is depicted as 
well. 
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Appendix 


Below we report on some values f(m, n ) and c(m, n ) that we have computed 

Note that 


by the algorithms for narrow strips described in Subsection 2.1 


Aichholzer’s results referred to in the captions of the tables have been obtained 
by a code that works for general point sets in the plane. 


n 

# unimodular triangulations of P 2 , n 

capacity 

2 

64 

1.500000 

3 

852 

1.622451 

4 

12170 

1.696380 

5 

182132 

1.747462 

6 

2801708 

1.784822 

7 

43936824 

1.813494 

8 

698607816 

1.836244 

9 

11224598424 

1.854774 

10 

181815529916 

1.870184 

11 

2964167665340 

1.883216 

12 

48580814410080 

1.894393 

13 

799696199314500 

1.904094 

14 

13212398835196240 

1.912597 

15 

218976668040908248 

1.920118 

16 

3639020246503687098 

1.926820 

17 

60616163842958990268 

1.932833 

18 

1011775545312594580868 

1.938260 

19 

16918718677672553292440 

1.943185 

20 

283368129709983000763876 

1.947675 

21 

4752924784523774226889308 

1.951787 

22 

79824154012907603962950312 

1.955568 

23 

1342199498257069824064033644 

1.959057 

24 

22592402326314503187343665228 

1.962288 

25 

380653341141186360494812030908 

1.965287 

375 

322379469795023850474295003353206721961082558051423084880244 

622493239840267086773982523966267988496124181049802588221265 

093372321162311900056540572387505992321869061059550930460251 

208699393384897025341791259853666381631934919090333735984559 

109034657337254608555310086159512051235323132000393390812675 

635858586268866021116515894234547255346166243854345190625384 

072567100646447104480055590853368584565574390095154916379988 

529322533797949278877762995003200075959780 

2.044130 


Table 3: Results for m = 2 (up to n — 15 by Aichholzer |2|).) 






n 

— 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

60 
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# unimodular triangulations of P, 3 jTI 

46456 

2822648 

182881520 

12244184472 

839660660268 

58591381296256 

4140106747178292 

295372308876234428 

21234538315776214604 

1535939689343151109944 

111655493479477379881272 

8150727077307189203809876 

597087996550303632801161860 

43871350204895836758556369212 

3231797978935266793268797809260 

238606105193380387756570932194588 

17651135152017098450035730535703808 

1308029292984065630362694842042395056 

97080539975603502667567153853690549804 

7215158047881650609075575773153609553148 

536905685776901371485436849505792415847140 

39997858254082097021224132017959794867440460 

2982752306685557862989393328648927558138800612 

222638546950211814977693932477091620801626551100 

16632293481947394846909242053460530217349594447732 

1243490745851056260557782821562507156418923407432920 

93034737749193459157244717739574844241159902101217660 

6965244882542454937020619818702059053741377290255068284 

521789556367416753405244328934612259884552790211207421244 

39111402471791798530057405675011481922023421912454706904712 


capacity 

1.722619 

1.785718 

1.829755 

1.861743 

1.886238 

1.905656 

1.921429 

1.934510 

1.945546 

1.954989 

1.963164 

1.970314 

1.976623 

1.982234 

1.987258 

1.991783 

1.995882 

1.999613 

2.003024 

2.006154 

2.009039 

2.011705 

2.014178 

2.016478 

2.018623 

2.020627 

2.022506 

2.024269 

2.025929 

2.027493 


140299226506605674595297010256056870425972772255987818065754 „ 0 r 12 o 6 
8248104578242160531202625665330141151767257855947224 


Table 4: Results for m = 3 (up to n = 10 by Aichholzer |2|.) 





Table 5: Results for m = 4 (up to n = 8 by Aichholzer f^].) 


4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


ft unimodular triangulations of P^ n capacity 

736983568 1.841066 

208902766788 1.880202 

61756221742966 1.908818 

18792896208387012 1.930751 

5831528022482629710 1.948080 

1835933384812941453312 1.962138 

584455230176565718869688 1.973785 

187686028049755013528577884 1.983601 

60685901262618326775192700244 1.991986 

19731268926382148037209063600412 1.999235 

6444884828545542240332780129017164 2.005567 

2113222804656668311309302902100087020 2.011147 

695163898467233943317499868644974218294 2.016103 

229316915701559537858641762255000442720116 2.020535 

75827610389461537709077484409103543784585710 2.024522 

25126215170054967918515172517611569017605311400 2.028130 

8341120564526486621411516194118239406742548614820 2.031409 

2773492141234111587868660829757660194821668552146720 2.034405 

923542836393484391634763396019705303948974557455679944 2.037151 

307927840785137540445620416718491375007439099052396877524 2.039680 

102788976776576952654837601573736966649525369087792041774754 2.042015 

34348007173983906887977536976143285540616055151406837464386952 2.044178 

11488735157368185493550332714289865763761071970295733672739011182 2.046188 

3846114798702837748429637799631405899720860204464245538618538080224 2.048061 

1288605098545745425026895981169173452565572617561190100930952999826626 2.049811 

432053057804721606144370507673814019049397228524243323674624721189513152 2.051449 

144960362040664009850310097114550102266012390064271163400409174321129109970 2.052986 

48666988587816100417429461522083834520238266167053837281579839713435133670880 2.054431 

16348321592766160525928861562545873132577437301251385744803248045014265044217096 2.055792 
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n 

# unimodular triangulations of P 5 , n 

capacity 

i 

252 

1.595455 

2 

182132 

1.747462 

3 

182881520 

1.829755 

4 

208902766788 

1.880202 

5 

260420548144996 

1.915513 

6 

341816489625522032 

1.941533 

7 

464476385680935656240 

1.961547 

8 

645855159466371391947660 

1.977388 

9 

913036902513499041820702784 

1.990240 

10 

1306520849733616781789190513820 

2.000871 

11 

1887591165891651253904039432371172 

2.009821 

12 

2747848427721241461905176361078147168 

2.017461 


Table 6: Results for m = 5 (up to n = 6 by Aichholzer J2|].) 


n 

# unimodular triangulations of Pe i7l 

capacity 

i 

924 

1.641958 

2 

2801708 

1.784822 

3 

12244184472 

1.861743 

4 

61756221742966 

1.908818 

5 

341816489625522032 

1.941533 

6 

1999206934751133055518 

1.965553 

7 

12169409954141988707186052 

1.984082 


Table 7: Results for m = 6. 








